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Abstract 


In this study, biodiesel fuel production from waste sunflower oil and viscosity optimization was carried out. During the pro- 
duction process, catalyst ratio, alcohol ratio and reaction temperature were determined as variable parameters. Transesferication 
method was used as the production method. During the production process, the use of NaOH catalyst and methyl alcohol was 
provided. Biodiesel production steps with the transesterification method were discussed in detail. A total of 27 different biodiesel 
fuels were obtained with a catalyst ratio varying between 0.03% and 0.07%, alcohol content between 15% and 25%, and reaction 
temperature between 50 ° C and 70 ° C. All biodiesel fuels were analyzed and their characteristics were determined. In the opti- 
mization process, catalyst ratio, temperature and alcohol ratio were considered as input parameters, and viscosity as output pa- 
rameters.Both 3D surface plots and 2D contour plots were developed using MINITAB 19 to predict optimum biodiesel viscosity. 
To predict biodiesel viscosity a quadratic model was created and it showed an R2 of 0.95 indicating satisfactory of the model. 
Minimum biodiesel viscosity of 4.37 was obtained at a temperature of 60, NaOH catalyst concentration of 0.07% and an alcohol 
ratio of 25%. At these reaction conditions, the predicted biodiesel viscosity was 4.247. These results demonstrate reliable predic- 


tion of the viscosity by Response surface methodology(RSM). 
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1. Introduction 


Due to the rapid increase in the world population and technologi- 
cal developments, the number of vehicles and energy need are in- 
creasing. Diesel engines are widely used in passenger and freight 
transportation, agricultural and industrial activities due to their high 
efficiency and durability [1]. These vehicles are dependent on fossil 
fuels to meet the energy needs [2,3]. But fossil energy sources are 
non-renewable energy sources. For this reason, it becomes difficult 
to provide it in a sustainable way and to ensure energy security [4]. 
Emissions from diesel engines bring negative effects along with it. 
This means that the exhaust emissions from vehicles will cause prob- 
lems for the environment and human health [5,6]. In order to close 
the gap between consumed energy and produced energy and to ob- 
tain sustainable clean energy, there has been a tendency towards al- 
ternative, renewable energy sources [7,8]. 
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An alternative fuel should be economically competitive with die- 
sel fuel, be environmentally safe and readily available. The most im- 
portant alternative fuels for diesel engines are biodiesel fuels under 
biomass energy [9,10]. Biodiesel fuel can be produced from vegeta- 
ble, animal or marine products. In addition, biodiesel can be pro- 
duced from waste oils [11,12]. Different methods are used for bio- 
diesel production. Biodiesel production can be realized by dilution, 
micro-emission, pyrolysis and transesterification methods [13,14]. 
The method of producing biodiesel as a result of the reaction of veg- 
etable and animal oils with alcohol is called transesterification. 
Transesterification method is widely used due to its low cost, mild 
reaction conditions, ease of production and properties close to stand- 
ard diesel fuel. This method also has disadvantages such as the dif- 
ficulty of separation processes, the risk of side reactions and the large 
amount of water waste [15,16]. There are almost no aromatic com- 
pounds, carcinogenic substances and sulfur in its structure. It can be 
used in standard diesel engine without any change. 
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Biodiesel fuel has many advantages. Biodiesel fuel has basic ad- 
vantages such as no risk of extinction, supporting agricultural devel- 
opment, reducing foreign dependency, and being environmentally 
friendly. The high cetane number decreases the knock value. Its high 
flash point increases transport and storage safety. It is a more lubri- 
cating fuel than standard diesel fuel and reduces wear values. Be- 
sides these advantages, there are also disadvantages. Its calorific 
value is lower than petrol based diesel fuel. This feature causes a 
slight decrease in power as a result of combustion in the engine, de- 
crease in engine performance and increase in specific fuel consump- 
tion. Its viscosity is higher than diesel fuel, it is affected more quickly 
by cold weather conditions. This situation is a factor limiting the use 
of biodiesel in engines in cold climates and in pure form 
[17,18,19,20]. 

In experimental studies, it is necessary to know the experiment 
design, parameters and what to expect from the test result in order to 
reach a correct result. Even if all these conditions are met, it may be 


necessary to make a large number of the same analysis or experiment. 


Reducing the number of experiments provides labor, time and cost 
savings. For this reason, a design determined with the correct param- 
eters and levels is a reason for preference [21]. Response Surface 
Method (RSM) is taking its place in current statistical methods 
[22,23]. RSM was developed by Box and Wilson in 1951 and was 
first applied to the chemical industry [24]. RSM is widely used in the 
formulation of a new product, improvement of existing product de- 
sign, process optimization, process development and improvement 
[25]. Response Surface Methodology can be used in experiments 
with at least two or more parameters. It is used for mathematical 
modeling and optimization of the relationship between outputs and 
inputs based on the results obtained from experiment combinations 
consisting of different levels of parameters. This method is also 
known as the experimental design method [21]. RSM uses experi- 
mental modeling techniques used to determine the relationship be- 
tween the system's response and the independent variables affecting 
it. It includes experimental strategies and optimization techniques to 
investigate the experimental space of process variables [26]. 

Anwar et al., have produced biodiesel from Australian native 
stone fruit by transesterification method. They took methanol oil 
molar ratio, catalyst ratio and temperature as variable parameters to 
use in the production process. Both 3D surface plots and 2D contour 
plots were developed using MINITAB 18 to estimate the optimum 
biodiesel yield. After optimization for maximum biodiesel yield, 
they found a methanol: oil molar ratio of 6: 1, KOH catalyst ratio 
0.5% and a reaction temperature of 55 C. They determined the bio- 
diesel yield produced under these reaction conditions as 95.9% [27]. 
Nayak et al., examined the optimization of methyl ester yield of bi- 
odiesel fuel produced from papaya oil using response surface 
method. They took temperature, catalyst amount, methanol / oil mo- 
lar ratio and reaction time as variable parameters. Based on the opti- 
mum condition, the predicted biodiesel yield was 99.9% and the ac- 
tual experimental value was 99.3% [28]. Latchubugata et al., opti- 
mized parameters such as temperature value, reaction time and 
methanol / oil molar ratio used in the process of biodiesel production 
from palm oil with the Response surface method. They produced bi- 
odiesel fuel at the parameter values reached as a result of optimiza- 
tion. As aresult of the analysis they applied to the produced biodiesel 
fuel, they found that the optimization reached the result with high 


accuracy [29]. The response surface method has been applied in dif- 
ferent working areas and successful results have been achieved. 
Giwvercin et al., using the Response Surface Method, optimized the 
cutting parameters that affect the surface roughness. They success- 
fully performed the determination of the parameter and the optimum 
value that most affected the results obtained in the experiments [26]. 
Ozmetin et al., conducted an experimental study on paint chemistry 
and optimized the results with the Response Surface Method. With 
the result they achieved by optimization, they reached the targeted 
values [21]. 

In this paper, biodiesel fuel production from waste sunflower oil 
and viscosity optimization were carried out depending on the alcohol 
ratio, catalyst ratio and temperature values. In the optimization per- 
formed with the response surface method, it has been concluded that 
biodiesel should be produced at 66.96 ° C temperature, 0.07% NaOH 
catalyst ratio and 25% alcohol ratio to obtain a minimum biodiesel 
viscosity of 4.227. 


2. Material and Method 


In this study, biodiesel fuel was produced from waste sunflower 
oil, depending on the alcohol ratio, catalyst ratio and temperature 
values. Transesterification method is used in the production of bio- 
diesel fuel. The biodiesel production process and the response sur- 
face method stages applied to the analysis results obtained were ex- 
amined in detail. 


2.1 Biodiesel production 


Transesterification method was used for the production of bio- 
diesel from waste sunflower oil. This method consists of 6 steps: 
mixing alcohol and catalyst, reaction, separation, alcohol removal, 
glycerin neutralization and methyl ester washing process. The stages 
of biodiesel production by the transesterification method can be seen 
schematically in figure 1 [30-34]. 


Waste sunflower oil 


Heating 
process 


Removal of 
alcohol 


Neutralization 
of glycerin 


Mixing 
alcohol and 
catalyst 


Washing 
methyl ester 


Figure 1. Biodiesel production process with transesterification 
method[30] 
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Beaker, glass bubble, cooler, density measuring device, separation characteristics of the materials used in biodiesel production are given 
funnel, heated magnetic stirrer with thermocouple, magnetic fish and __ in table 1. 
precision scales were used for biodiesel production. The technical 


Table 1. Technical properties of materials used in biodiesel production 
Material name Waste sunflower oil Sodium hydroxide (NaOH) Methyl alcohol (CH30H) 


Density (g/cm?) 0.922 2.13 0.790-0.793 
Molecular weight (g/mol) - 40 32.04 
Boiling point - - 64-65 
Melting point - 319-322°C - 
Resolution 1090 g/l = = 
Refractive index - - 1.328-1.331 


Table 2. Technical properties of viscosity device 
Device Name Brand / Model Measuring Range Temperature Range Sample Volume 
Viscosity Omnitek D445 0.15-25,000 mm?/s @ 40 °C 15-150 °C 8-16 ml 


Waste sunflower oil was heated up to 80 ° C and then filtered. In 
order to evaporate the water molecules in the oil, it was kept at 120 ° 
C for about 1 hour. A homogeneous solution was prepared with 
NaOH catalyst and methyl alcohol and mixed with waste sunflower 
oil. Twenty-seven different biodiesel production was carried out, 
with catalyst ratio varying between 0.03% and 0.07%, alcohol con- 
tent between 15% and 25%, and reaction temperature between 50 ° 
Cand 70 ° C. The reaction time was kept constant at 1.5 hours in all 
production processes. The setup where the reaction taking place dur- 
ing the production process takes place is shown in Figure 2. 


Figure 2. Reaction setup 


After the reaction, the separation of glycerin from biodiesel took 
approximately 10 hours. The glycerin separation process in the bio- 
diesel production process is shown in Figure 3.The fatty acids, cat- 
alysts and unreacted alcohols contained in the biodiesel obtained as 
a result of the reaction were removed by washing. The washing 
process was repeated until the fuel cleared. For this study, 5 repeti- 
tive washing processes were carried out. The washing process is 
shown in figure 4. After the washing process, drying and filtering Figure 4. Biodiesel fuel washing process 
process was applied to biodiesel fuel again. 
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2.2 Design of Experiments 


Response surface methodology was used for statistica analysis of 
the experimental data using the MINITAB 19 software. The Box- 
Behnken is one of the most commonly used response surface meth- 
odology designs. This design was used for statistical analaysis and 
designing of this experiment. The Box-Behnken design matrix was 
utilised to find the optimum conditions for minimum biodiesel vis- 
cosity. The experimental optimization was achieved via ANOVA 
(analaysis of variance) using MINITAB 19 software. The effect of 
process factors such as alcohol, NaOH catalyst concentration, and 
temperature were tested. Using these three factors at three levels re- 
quired a total of 15 runs for identifying the optimum conditions for 
biodiesel viscosity. The ranges, coded symbols and levels of the fac- 
tors are shown in Table 3. The design matrix fort he three factors 
was varied at three levels, namely -1, 0 and +1. 


Table 3. Range and levels coded for independent factors 


Alcohol ratio ranged from 15% to 25%, NaOH catalyst concen- 
trations were 0.03 — 0.07% and the reaction temperature was varied 
from 50 °C to 70 °C. The response factor (biodiesel viscosity) was 
correlated to the parameters using a full quadratic model. The gen- 
eral form of full quadratic model is expressed as follows, 


Y = Bo + By X, + BoX2 + B3X3 + By 2X1X2 + By 3X1 X3 + 


By 3X2X3 + By aX,” + By 2X,” + B3 3X3" (1) 


where Y is the predicted biodiesel viscosity; Bp is a constant; B,, 
B,, and Bz are regression coefficients; B,,, By,2, By3, and Bz, 
are quadratic coefficient; and X,, Xz, and X3 are independet vari- 
ables. 


3. Simulation Results 


The results of the Box-Behnken design model to optimize bio- 


Range and Levels 


diesel viscosity parameters are shown in Table 4. In the experimental 


Factors Unit phe results, the viscosity of biodiesel ranged from 4.376 to 5.478 mm2/s. 
-1 0 +1 This design matrix also show the run order, experimental viscosity 
values and predicted viscosity values. To avoid systematic errors, all 
Alcohol (%) A 15 20 25 run orders were randomised. 
Nach cays des Cc 0.03 0.05 0.07 
concentration 
Temperature °C T 50 60 70 
Table 4. Experimental matrix and Box-Behnken results 
Exp. Run n : Temp. Alcohol ea Biodiesel Viscosity 
Number Order (°C) (%) (%) Experimental Predicted 
1 5 -1 0 -1 50 15 0,05 5,058 5,00575 
2 7 -1 0 1 50 25 0,05 4,518 4,53250 
3 10 0 1 -1 60 15 0,07 4,627 4,56562 
4 8 1 0 1 70 25 0,05 4,442 4,49425 
5 6 1 0 -1 70 15 0,05 5,090 5,07550 
6 4 1 1 0 70 20 0,07 4,423 4,49887 
7 14 0 0 0 60 20 0,05 4,816 4,81600 
8 11 0 -1 1 60 25 0,03 4,759 4,82037 
9 15 0 0 0 60 20 0,05 4,816 4,81600 
10 2 1 -1 0 70 20 0,03 5,478 5,36437 
11 3 -1 1 0 50 20 0,07 4,453 4,56662 
12 9 0 -1 -1 60 15 0,03 5,429 5,55713 
13 12 0 1 1 60 25 0,07 4,376 4,24787 
14 13 0 0 0 60 20 0,05 4,816 4,81600 
15 1 -1 -1 0 50 20 0,03 5,341 5,265.12 


The predicted biodiesel viscosity values were obtained from 
Minitab software version 19.0 using a quadratic regression model by 
means of response surface methodology (RSM) analysis of experi- 
mental data. Minitab 19 program was used to calculate each param- 
eter and the effects of their interactions with other parameters. Bio- 
diesel viscosity was correlated with other parameters using the quad- 
ratic regression model shown in Equation (2). 
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V = 7.22 — 0.0302T — 44C + 0.059A + 0.000435T2 + 
161C? — 0.0033A? — 0.209TC — 0.00054TA + 1.047KA (2) 


Here, V is response, C is catalyst concentration, T is a reaction 
temperature, and A represents alcohol ratio. 
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Viscosity 
w 


—— Experimental ——Predicted 


1 2 3 4 5 6 7 8 9 10 


Number of Experiment 


5. Experimental and RSM fitting 


Comparison of experimental and predicted biodiesel are shown in 
Fig. 5. It is seen that there is a sufficient correlation between RSM 
predictive values and experimental values confirming the accepta- 
bility of the model. 


(response is Viscosity; « = 0,05) 


2,571 
T 


Factor Name 
Temperature 
B Catalyst 
\c Alcohol 


i) 1 2 : 3 4 5 6 rd 8 9 
Standardized Effect 


Figure 6. Effects of factors 


Table 5. Regression coefficients 


Standard 


Term. Coefficients T-Value P-Value 
Errors 

Constant 4,8160 0,0746 64,52 0,000 
T 0,0079 0,0457 0,17 0,870 
C -0,3910 0,0457 -8,55 0,000 
A -0,2636 0,0457 -5,77 0,002 
TXT 0,0435 0,0673 0,65 0,546 
CXC 0,0643 0,0673 0,95 0,383 
AXA -0,0825 0,0673 -1,23 0,275 
TXC -0,0417 0,0646 -0,65 0,547 
TXA -0,0270 0,0646 -0,42 0,694 
KXxA 0,1048 0,0646 1,62 0,166 


In order to investigate the effects on biodiesel viscosity, linear, 
quadratic and interaction effects of parameters were taken into ac- 
count. Table 5 and Figure 6 show the importance of these parame- 
ters in terms of the probability value (p-value). It also shows the 
obtained regression coefficients and calculated T-values. In the 
model, positive coefficients T, T2, K and KA had a positive effect 
on biodiesel viscosity, while A, A2, TC, TA and C had negative ef- 
fects on biodiesel viscosity. Analysis of variance (ANOVA) was 
used to determine the importance and appropriateness of the quad- 


ratic model. 


Table 6. ANOVA results for biodiesel viscosity 


Degree of F Sum of Squ 


Mean squar 


Source F-Value P-Value Remarks 
reedom ares e 
Model 9 1,88910 0,20990 11,22 0,008 Significant 
Linear 3 1,77953 0,59318 31,71 0,001 Significant 
T-Temperature 1 0,00050 0,00050 0,03 0,877 Not significant 
C-Catalyst 1 1,22305 1,22305 65,37 0,000 Highly significant 
A-Alcohol 1 0,55599 0,55599 29,72 0,003 Significant 
Square 3 0,05580 0,01860 0,99 0,467 Not significant 
T? 1 0,01271 0,01271 0,68 0,447 Not significant 
Cc 1 0,02329 0,02329 1,24 0,315 Not significant 
A? 1 0,01674 0,01674 0,89 0,388 Not significant 
TC 3 0,05378 0,01793 0,96 0,480 Not significant 
TA 1 0,00697 0,00697 0,37 0,568 Not significant 
CA 1 0,00292 0,00292 0,16 0,709 Not significant 
Lack-of-Fit 1 0,04389 0,04389 2,35 0,186 Not significant 
Pure Error 5 0,09354 0,01871 
Total 3 0,08358 0,02786 5,59 0,155 
R2=0.9575 2 0,00996 0,00498 
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Table 6 shows the significance of the individual terms and their in- 
teraction on the selected response. The P value is used to check the 
significance of each regression coefficient by representing the error 
probability. The interaction effect of each cross product can be re- 
vealed through the p value[27]. It is found, C (Catalyst concentra- 
tion), A (Alcohol) have significant effects on biodiesel viscosity. It 
is seen that C value has the lowest p value (0.000) and the highest F 
value (73.16) according to all other parameters. These results show 
that the C value is the most important parameter in biodiesel viscos- 
ity. According to the regression model in Equation (1), A has a pos- 
itive effect and both C and reaction temperature (T) have negative 
effects on biodiesel viscosity. This implies that increasing A will in- 
crease the viscosity of the biodiesel. However, increase in C and T 
will decrease the viscosity of biodiesel. The ANOVA results showed 
that the linear term of T with p value was not significant (more than 
0.05) and its quadratic term T2 with p value also was significant 
(more than 0.05). R2 also shows good correlation between independ- 
ent parameters. In this study, R2 was found to be 95.75% and the 
corrected coefficient of determination (Adj. R2) was found to be 
88.11%. This means that the model explains 95.75% of the variation 
in experimental data. As a result, the regression model developed for 
biodiesel viscosity was valid and showed a satisfactory experimental 
relationship between response and parameters. 


Viscosity 5,0 


25 


20 Alcohol 


0,075 


a) 


Alcohol 


3.1 Interaction effect of Alcohol and Temperature 


The interaction effect of alcohol A and catalyst concentration, C 
on biodiesel viscosity in both the 3D surface plot and the contour 
plot are shown in Figure 3. With an increase of catalyst concentration 
0.07(highest) and alcohol 25% (highest) biodiesel viscosity decrease. 
The minimum viscosity of biodiesel value of 4.227 mm2/s was 
found for NaOH 0.07 % (Run 12). Table 3 design matrix indicated 
that lowest NaOH concentration at 0.03% and mid-value of alcohol 
ratio at 20% resulted in highest biodiesel viscosity. When the Alco- 
hol ratio remains unchanged at 25% and catalyst concentration is at 
lowest value of 0.03, the biodiesel viscosity decreases 4.759 mm2/s 
(Run 11). When the alcohol ratio was reduced to 15 % (lowest level), 
and with the highest value of catalyst concentration of 0.07%, the 
biodiesel viscosity was found to be 4.627 (Run10). Again, at alcohol 
ratio of 25%, and with the mid-value of catalyst concentration of 
0.05%, the viscosity was found to be 4.442% mm2/s (Run 7). On the 
other hand, when the alcohol ratio was reduced to 15, and with the 
mid-value of catalyst concentration of 0.05% the yield rose up to 
5.09 mm2/s (Run 6). Alcohol ratio affected total biodiesel viscosity. 
ANOVA from Table 6 confirmed that both A and C interaction were 
significant. The 2D contour plot with A and C interaction along with 
biodiesel viscosity is shown in Figure 7. It is easy to identify the op- 
timum operating conditions and the related response values (viscos- 
ity) through the 2D contour plot. Therefore, both A and C are signif- 
icant for lower biodiesel viscosity. 


Viscosity 
< 44 


a Tt | 
LO PEe 
NOOO 


virtues 
co) 
p 


0,05 


Catalyzer 
b) 


Figure 7. Interaction effect of Alcohol and catalyst concentration on bioediesel viscosity a) 3D surface plot, b) Contour plot 


4. Conclusions 


In biodiesel production process, catalyst ratio, alcohol ratio and 
reaction temperature were determined as variable parameters. Trans- 
esferication method was used as the production method. During the 
production process, the use of NaOH catalyst and methyl alcohol 
was provided. Biodiesel production steps with the transesterification 
method were discussed in detail. A total of 27 different biodiesel 
fuels were obtained with a catalyst ratio varying between 0.03% and 
0.07%, alcohol content between 15% and 25%, and reaction temper- 
ature between 50 ° C and 70 ° C. All biodiesel fuels were analyzed 
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and their characteristics were determined. A response surface meth- 
odology based on Box—Behnken design matrix was applied to 
achieve the optimum biodiesel viscosity. Three main parameters 
were changed separately at different intervals to estimate the bio- 
diesel viscosity in this matrix. Based on the results, optimum bio- 
diesel viscosity were found to be Alcohol ratio of 25%, catalyst con- 
centration 0.07%, and a reaction temperature of 66.96 °C. The min- 
imum biodiesel viscosity under such conditions was 4.376 mm2/s, 
which also confirmed the RSM model prediction of 4.227 mm2/s. 
ANOVA statistics of this study confirmed that catalyst concentration 
ratio has the most significant effect on the biodiesel viscosity, 
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whereas reaction temperature does not seem to have any significant 
impact. 


Nomenclature 


ANOVA 
RSM 


Analysi of Variance 
Response Surface Method 
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